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Type I planetary nebulae (PNe) have high He/H and N/0 ratios and are 
thought to be descendants of stars with initial masses of ~ 3 - 8M 0 . These 
characteristics indicate that the progenitor stars experienced proton-capture nu- 
cleosynthesis at the base of the convective envelope, in addition to the slow 
neutron capture process operating in the He-shcll (the s-process). We compare 
the predicted abundances of elements up to Sr from models of intermediate-mass 
asymptotic giant branch (AGB) stars to measured abundances in Type 1 PNe. 
In particular, we compare predictions and observations for the light trans-iron 
elements Se and Kr, in order to constrain convective mixing and the s-process 
in these stars. A partial mixing zone is included in selected models to explore 
the effect of a 13 C pocket on the s-process yields. The solar-metallicity mod- 
els produce enrichments of [(Se, Kr)/Fe] < 0.6, consistent with Galactic Type 
I PNe where the observed enhancements are typically < 0.3 dex, while lower 
metallicity models predict larger enrichments of C, N, Se, and Kr. O destruction 
occurs in the most massive models but it is not efficient enough to account for 
the > 0.3 dex O depletions observed in some Type I PNe. It is not possible to 
reach firm conclusions regarding the neutron source operating in massive AGB 
stars from Se and Kr abundances in Type I PNe; abundances for more s-process 
elements may help to distinguish between the two neutron sources. We predict 
that only the most massive (M > 5M 0 ) models would evolve into Type I PNe, 
indicating that extra-mixing processes are active in lower-mass stars (3-4M 0 ), if 
these stars are to evolve into Type I PNe. 

Subject headings: stars: AGB and post-AGB stars — planetary nebulae: general 
— nuclear reactions, nucleosynthesis, abundances 


1. Introduction 

After the thermally-pulsing AGB (TP-AGB) phase is terminated, low to intermediate 
mass stars (~ 0.8 to 8M 0 ) evolve to become post-AGB stars and possibly planetary nebulae 
(PNe), if the ejected envelopes have sufficient time to become ionized by the hot central 
stars before dissipating into the interstellar medium. For recent reviews of TP-AGB and 
post-AGB stars see Herwig (2005) and van Winckel (2003), respectively. The illuminated 
PN is comprised of material from the deep convective envelope, hence nebular abundances 
should reveal information about the efficiency of mixing events and chemical processing that 
took place during previous evolutionary phases, in addition to the initial composition of the 
parent star (Dopita et al. 1997; Karakas & Lattanzio 2003a). 



Briefly, during the TP-AGB phase the He-burning shell becomes thermally unstable 
every 10 4 years or so, depending on the core mass. The energy from the thermal pulse 
(TP) drives a convective zone in the He-rich intershell, that mixes the products of He- 
nucleosynthesis within this region. The energy provided by the TP expands the whole star, 
pushing the H-shell out to cooler regions where it is almost extinguished, and subsequently 
allowing the convective envelope to move inwards (in mass) to regions previously mixed 
by the flash-driven convective zone. This inward movement of the convective envelope is 
known as the third dredge- up (TDU), and is responsible for enriching the surface in 12 C 
and other products of He-burning, as well as heavy elements produced by the s process. 
Following the TDU, the star contracts and the H-shell is re-ignited, providing most of the 
surface luminosity for the next interpulse period. In intermediate-mass AGB stars with 
initial masses > 4 M 0 , the base of the convective envelope can dip into the top of the H- 
burning shell, causing proton-capture nucleosynthesis to occur there (hot bottom burning, 
HBB; see Lattanzio et al. 1996, for a review). 

Abundances can be obtained from PN spectra for a number of elements including C, N, 
O, S, Cl and the noble gases He, Ne, and Ar (Allcr & Czyzak 1983; Kingsburgh & Barlow 
1994; Dopita et al. 1997; Stanghellini et al. 2000; Leisy & Dennefeld 2006); and heavy 
elements synthesized by the s process, including Ge, Se, Kr, Xe, Ba, and Br (Pequignot & 
Baluteau 1994; Dinerstein 2001; Sterling et al. 2002; Sharpee et al. 2007; Sterling et al. 2007; 
Sterling & Dinerstein 2008, hereafter SD08). Whereas the metallicity of a star is usually 
taken from the abundance of iron, for ionized nebulae oxygen is generally the primary index 
of metallicity. However, O is not an ideal metallicity index for PNe, because it can be 
destroyed by CNO cycling and/or synthesized during helium burning. Consequently, Ar is 
sometimes used in place of O (Leisy & Dennefeld 2006, SD08), under the assumption that it 
remains unchanged by AGB nucleosynthesis. Zinc is a better indicator of iron-group element 
abundances in nebulae than Fe itself, because Zn does not condense into dust as easily as 
Fe (York & Jura 1982; Sembach et al. 1995; Welty et al. 1999) and can be observed in 
PNe (Dinerstein & Geballe 2001; Dinerstein 2004). However, Zn is at the beginning of the 
s-process chain and might be produced to some extent in AGB stars. In order to interpret 
the elemental abundances in PNe it is important to determine which elements are produced, 
destroyed or left unaltered by AGB nucleosynthesis. 

The set of elements heavier than Fe that can be detected in PNe include all even- 
numbered elements from Z = 30 - 36, which lie at the light-element end of the s-process 
distribution, just below the first s-process peak at Sr, Y, and Zr (Z — 38 - 40). In Karakas 
et al. (2007) we attempted to use Ge abundances derived from a small sample of PNe (Sterling 
et al. 2002) to try and constrain the amount of TDLI mixing after the final few TPs. There a 
hint was found that efficient TDLI at the tip of the AGB may be required in order to match the 
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abundances; however, stronger conclusions were precluded by model uncertainties, the small 
PNe sample size, and the large uncertainties of the Ge abundance determinations. In the 
present study we proceed to analyze Se and Kr abundances from a much larger sample of 120 
Galactic PNe from SD08, where the .s- pro cess enrichments were investigated for correlations 
with PN morphology and other nebular and stellar characteristics, as outlined in detail in 
§2. We also use PNe abundances of C, N, and O along with the abundances of post-AGB 
stars for comparison to the stellar models. 

The composition of Type I PNe provides important constraints on nucleosynthesis and 
mixing in intermediate-mass AGB stars. These PNe tend to exhibit bipolar morphologies and 
have high He/H and N/O ratios characteristic of proton-capture nucleosynthesis (Stanghellini 
et al. 2006), suggesting they are descendants of AGB stars with initial masses ~3 - 8 M & 
(Peimbert 1978; Kingsbnrgh & Barlow 1994; Stanghellini et al. 2006). Low-mass AGB 
stars, on the other hand, are defined as having initial masses ~1 - 3 M 0 and tend to show 
carbon and s-process element enrichments (Busso et al. 2001; Travaglio et al. 2004). In the 
Milky Way, the spatial distribution and kinematics of Type I PNe indicate that they are a 
young population (e.g., Corradi & Schwarz 1995). This supports their identification with 
intermediate-mass stars, which have relatively short evolutionary lifetimes. In the Galaxy, 
intermediate-mass AGB stars are difficult to identify owing to a lack of reliable distances. 
This has resulted in a paucity of observational evidence for constraining stellar models, 
and especially the efficiency of the TDLT in intermediate-mass AGB stars, which is still not 
well determined (Karakas et al. 2002; Ventura et al. 2002; Herwig 2004; Stancliffe et al. 
2004). Garcfa-Hernandez et al. (2006) identified several Galactic intermediate-mass AGB 
stars within a sample of OH/IR stars (i.e., bright O-rich giants with large infrared excesses). 
The large enhancements of Rb found by these authors, combined with the fact that these 
stars are O-rich, support the prediction that HBB and an efficient TDU has occurred in 
these stars. At lower metallicities, constraints are provided by observations of luminous 
O-rich AGB stars in the Large and Small Magellanic Clouds (LMC and SMC, respectively) 
that are rich in Li and s-process elements (Wood et al. 1983; Smith & Lambert 1989; Plez 
et al. 1993). 

We begin in §2 with a summary of observational results for n-capture elements in PNe. 
We review the numerical method and present the stellar models in §3, and outline the results 
in §4. In §5, we discuss our findings and their implications, and in §6 we summarize our 
conclusions. 
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2. Observational Constraints 

The low cosmic abundances of n-captnre elements cause their spectroscopic features to 
be very weak. For that reason, trans-iron elements were not identified in the spectrum of a 
PN until 1994 (Peqnignot & Baluteau 1994). However, in the last few years there has been 
considerable progress in the detection of emission lines from and derivation of abundances 
for trans-iron elements in PNe (e.g., Dinerstein & Geballe 2001; Sterling et al. 2007; Sharpee 
et al. 2007, SD08). 

The most comprehensive study is that of SD08, who conducted the first large-scale 
survey of n-capture elements in PNe. They determined Kr and Se abundances in 120 PNe 
from the near-infrared emission lines [Kr III] 2.199 /mi and [Se IV] 2.287 /mi identified by 
Dinerstein (2001), and corrected for the abundances of unobserved Se and Kr ions using 
analytical formulae derived from a grid of photo-ionization models (Sterling et al. 2007). 
This survey increased the number of PNe with determined n-captnre element abundances 
by nearly an order of magnitude, enabling a detailed study of s-process enrichments in PNe 
as a population, as well as an assessment of enrichment patterns in different classes of PNe. 
In order to normalize the abundances of the s-process products to the initial metallicity of 
each PN, the Se and Kr abundances were compared to those of O and Ar reported in the 
literature. Oxygen was utilized as a reference element for PNe with low-mass progenitors, 
since its abundance is generally the most reliably determined in PNe, while Ar was used 
for objects with higher-mass central stars that may have experienced O destruction during 
HBB. Overall, 41 of the 94 PNe with derived Se and/or Kr abundances or meaningful upper 
limits were found to be enriched in these elements, with the average [Se/(0, Ar)]= 0.31 and 
[Kr/(0, Ar)]= 0.98 1 . SD08 interpreted these enrichments as evidence for in situ s-process 
nucleosynthesis and dredge-up in PN progenitor stars. 

SD08 observed significant numbers of PNe with intermediate-mass progenitor stars, as 
evidenced by their Type I compositions and/or bipolar morphologies (29 and 28 objects 
respectively; see their Table 15). They found that Type I and bipolar PNe exhibit smaller 
Se and Kr enrichments on average than other PNe. In Figure 1 we display [Se/(0, Ar)] and 
[Kr/(0, Ar)], separated into 0.1 dex bins, for non-Type I, Type I, and bipolar PNe. The 
distribution of Se and Kr enrichments is clearly skewed to smaller values for Type I PNe 
than other objects, and to a lesser extent the same is true of bipolar PNe. Notably, all of 
the Type I PNe observed by SD08 exhibit modest (less than a factor of two, or 0.3 dex) or 
no Se and Kr enrichments. These values are significantly lower than the average values for 


1 We use the notation [Se/(0, Ar)] and [Kr/(0, Ar)] of SD08 to emphasize the use of different reference 
elements for PNe with different progenitor masses, where [X/Y] = logio(X/Y) — logio(X/Y) 0 . 
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non-Type I PNe: 0.36 dex for [Se/(0, Ar)], and 1.02 dex for [Kr/(0, Ar)]. In contrast, while 
some bipolar PNe are among the objects with the lowest [Se/(0, Ar)] and [Kr/(0, Ar)], 
others display enrichments approaching a factor of 10. 

Figure 2 shows [Se/(0, Ar)] and [Kr/(0, Ar)] as a function of the logarithmic He/H, 
N/O, and C/O abundances of PNe from the SD08 sample, including upper limits (open 
symbols) for objects without detected Se and/or Kr emission. This figure not only illustrates 
the lower average Se and Kr abundances of Type I PNe relative to other objects, but also that 
upper limits to the Se and Kr abundances in many Type I PNe allow for at most marginal 
s-process enrichments. SD08 found that s-process enrichments are positively correlated with 
C/O in non-Type I PNe (bottom panels of Figure 2), a result that agrees with studies of AGB 
(Smith & Lambert 1990; Abia et al. 2002) and post- AGB stars (van Winckel 2003). Type I 
PNe often exhibit low C/O ratios (< 1), in agreement with predictions from intermediate- 
mass stars with HBB. HBB converts dredged-up 12 C into 14 N, thereby preventing or delaying 
the formation of C-rich AGB stars (Frost et al. 1998). On the other hand, some bipolar 
PNe display significant enrichments of both C and s-process nuclei, as well as non-Type I 
compositions. These objects may be descendants of low-mass binary star systems in which 
the bipolar morphology is a result of binary interactions (Soker 1997; Balick & Frank 2002). 
This suggests that some bipolar PNe are not descendants of intermediate-mass stars, and 
that morphology is not as reliable an indicator of progenitor mass as is chemical composition. 

Sharpee et al. (2007) derived n-capture element abundances in five PNe from high 
resolution optical spectra. They identified emission lines of several n-capture elements, 
including Br, Kr, Xe, Rb, Ba, and Pb in each of these objects, and derived abundances for 
Br, Kr, and Xe. They found that Kr and Xe are enriched in three of the five observed PNe, by 
0.3 to 0.9 dex. The two objects in their sample which do not exhibit s-process enrichments, 
NGC 2440 and IC 2501, are both Type I PNe according to the classification scheme of 
Peimbert (1978). The [Kr/Ar] and [Xe/Ar] values of these two objects are consistent within 
the uncertainties with the solar values, except for Xe in NGC 2440, which may be mildly 
subsolar. These results for [Kr/Ar] are consistent with the findings of SD08 for their larger 
sample of Type I PNe. 

The near-infrared [Kr III] and [Se IV] lines utilized by SD08 have also been detected 
in two PNe belonging to metal-poor Milky Way dwarf satellite galaxies (Wood et al. 2006, 
Dinerstein et ah, in preparation). Such observations can provide important constraints 
on AGB models and nucleosynthetic yields for subsolar metallicities. One of the observed 
objects, LMC SMP 62, is N-rich (N/O ~ 0.5) and hence considered a Type I PN in the 
context of the LMC, where the threshold value of N/O for Type I objects is taken to be 
lower than for Galactic PNe (Dopita & Meatheringham 1991; Leisy & Denncfeld 2006). 
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While the [Kr III] line is not detected (and does not provide a useful abundance limit), 
the [Se IV] line yields a Se abundance ~ 15% solar (Dinerstein et ah, in preparation). Not 
only is Se not enriched, it is actually deficient relative to the ce-elements; the abundances of 
O, Ne, S, and Ar in SMP 62 are ~ 30-40% solar (Leisy & Dennefcld 2006; Bernard-Salas 
et al. 2008) 2 . This is reminiscent of the abundance pattern in LMC red giants, where light 
s-process elements such as Y and Zr are markedly deficient relative to iron and some a 
elements (Pompeia et ah 2006; Mucciarelli et ah 2008). In contrast, Hen 2-436, a PN in the 
Sagittarius dSph galaxy, displays Se and Kr abundances as high as those for some of the 
most enriched Galactic objects shown in Figure 1. However, Hen 2-436 is C-rich rather than 
N-rich, and therefore is a non-Type I PN, with an estimated initial stellar mass of ~ 1.2 M© 
(Dudziak et ah 2000; Zijlstra et ah 2006). 

Zinc is detectable in PNe via the [Zn IV] 3.625 //rn fine-structure line identified by 
Dinerstein & Geballe (2001). This line has recently been observed in about a dozen PNe 
belonging to the Milky Way and its dwarf satellites (Dinerstein et ah 2007, Dinerstein, 
Geballe, & Sterling, in preparation, hereafter DGS08). Two of these objects exhibit Type I 
compositions: LMC SMP 62, discussed above; and the Galactic PN M 1-40 (Gorny et ah 
2004). The gaseous Zn abundance in M 1-40 is ~ 80% solar, which is consistent with its 
approximately solar Ne and S abundances. Even allowing for modest depletion into dust, 
which would raise the total Zn abundance by 0.1-0. 2 dex (assuming conditions typical of 
the warm ISM; see §1), it does not appear that Zn is significantly enriched in M 1-40. In 
LMC SMP 62, the measured Zn abundance is lower relative to solar than the abundances of a 
elements such as O, Ne, and S (Dinerstein et ah 2007) by a factor very similar to the deficiency 
in Se. In this context, it is interesting that several of the iron-group elements, including Ni 
and Cu (which are adjacent to Zn), display behavior similar to light n-capture elements, 
which are also deficient relative to Fe in LMC giants (Pompeia et ah 2006; Mucciarelli et ah 
2008). 


3. The stellar models 

The numerical method we use has previously been described in detail (Lugaro et ah 
2004; Karakas et ah 2006). We first compute the stellar structure using the Mt Stromlo 
Stellar Structure code (Lattanzio 1986), where each model is evolved from the zero-age main 


2 Using Spitzer spectra, Bernard-Salas et al. (2008) find higher Ne/H than Leisy & Dennefeld (2006). 
However, they adopt a higher value for the solar Ne abundance, and consequently their ratio to solar is 
similar. They also derive a low value for S/H, but the significance of this result is unclear, since these 
authors also find that Galactic PNe have subsolar sulfur abundances (Pottasch & Bernard-Salas 2006). 



sequence to near the tip of the TP-AGB. For the models considered in this study, we do not 
evolve to the post-AGB phase due to convergence difficulties that occur when the envelope 
mass is reduced below ~ 1.5M 0 ; see Karakas & Lattanzio (2007) for more details. The 
structure is then used as input into a post-processing nucleosynthesis code where we obtain 
abundances for many more species (up to 166) than are included in the structure model (6 
species). 

The stellar models we computed (see Table 1) cover a mass range of 3.O-6.5M 0 and 
metallicities 0. 2-1.0 Z & . The metallicities were chosen to reflect the composition of most 
of the PNe in the SD08 sample, where O and Ar abundances were used as metallicity 
indicators. We include a 2.5 M 0 , Z = 0.008 model as an example of a low-mass, low- 
metallicity progenitor that would lead to a non-Type 1 PN composition. All models use 
scaled solar abundances from Anders & Grevesse (1989) except for the Z = 0.012 models, 
which are computed with the revised solar elemental abundances from Asplund, Grevesse, & 
Sauval (2005) for comparison. Reimers (1975) mass loss was used on the first giant branch, 
with the parameter // = 0.4, and the Vassiliadis & Wood (1993) mass loss on the AGB. For 
comparison, we include results for a 5 M 0 , Z = 0.02 model computed with Reimers mass 
loss on the AGB, with the parameter r/ = 3.5. 

In Table 1 we summarize structural information about the models. Most of the structure 
calculations are taken from Karakas & Lattanzio (2007), with the Z = 0.012 models discussed 
previously in Karakas et al. (2007). In Table 1 we present the initial mass and metallicity, 
Z ; the total number of TPs computed; the maximum temperature in the He-shell, Tne™) the 
maximum temperature at the base of the convective envelope, T^g X ; the total mass dredged 
into the envelope during the TP-AGB, M dred ; the final envelope mass M env ; and whether 
or not HBB occurs. In Table 1, masses are given in solar units, temperatures in millions of 
kelvins, metallicity Z is the fraction by mass of metals, and elemental and isotopic ratios are 
in terms of number ratios. 


3.1. Type I planetary nebulae and HBB 

Observations suggest that the Type I PN progenitor mass range is ~ 3 to 8 M 0 (e.g., 
Peimbert 1990; Stanghellini et al. 2006). For the purpose of this study we consider models 
with masses between 3 to 6.5 M 0 , noting that more massive AGB stars would likely evolve 
too quickly to form observable PNe. This mass range encompasses the full range of AGB 
nucleosynthetic behavior: (1) models that experience the TDU but not HBB, and (2) models 
that undergo the second dredge-up (SDLI) and HBB in addition to the TDLL From inspection 
of Table 1, the 2.5, 3, and 4 M 0 models fall into the first category while the 5, 6, and 6.5 M 0 
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models fall into the latter category. The more massive objects that experience both the SDU 
and HBB are often defined as intermediate- mass AGB stars (see Karakas et al. 2006; Herwig 
2005). 

Only models with HBB (i.e., models with initial masses m > 5 M 0 ) would evolve into 
Type I PNe. Hence efficient extra mixing between the base of the convective envelope and 
a region hot enough to allow some ON cycling is presumably operating to produce the 
observed Type I nucleosynthethic signature of high He/H and N/O ratios in the lower-mass 
~ 3-4 M 0 stars. The physical mechanism that could cause such efficient extra mixing is 
not known, although stellar rotation is one possible candidate (see below). It is probably 
not thermohaline mixing which is only efficient in stars less massive than 2 M 0 (Eggleton 
et al. 2006; Charbonncl & Zahn 2007). The idea that extra-mixing processes are operating 
in low-mass carbon-rich AGB stars is supported by observational evidence, including lower 
than predicted 12 C/ 13 C ratios observed in N-type AGB stars (Abia & Isern 1997). The C, 
O, and Al isotopic ratios found in meteoritic silicon carbide (SiC) and oxide grains from 
AGB stars also point toward some sort of extra-mixing process operating in the parent stars 
(Busso et al. 1999; Nollett et al. 2003; Zinner et al. 2006; Zinner 2008). In these cases, the 
mixing should not be so efficient to prevent the formation of a carbon-rich atmosphere. 

Some fraction of Type I and bipolar PNe presumably formed from low-mass stars via 
binary evolution (Soker 1997; Balick & Frank 2002; Moe & De Marco 2006), although it is 
unclear how this would lead to the high He/H and N/O ratios observed. Two-dimensional 
simulations of a star with a point-mass companion in a very close orbit suggested that 
close binary evolution has little effect on the deep interior structure, despite producing large 
distortions to the outer layers (Deupree & Karakas 2005). However, these simulations did 
not go beyond the point where the stars begin to fill their Roche lobes. Rapid stellar 
rotation in combination with magnetic fields has also been proposed as a mechanism to 
shape bipolar PNe and to produce the overabundances of He and N in single intermediate- 
mass progenitors (M > 3 M 0 ) (Calvet & Peimbert 1983; Gorny et al. 1997; Garcia-Segura 
et al. 1999; Dobrincic et al. 2008). Further modeling is required to address how effectively 
binary interactions, rotation, and/or magnetic fields can affect the internal structure and 
nucleosynthesis of AGB stars. 

In Fig. 3 we plot the temporal evolution of the core, and the temperature in the He- 
shell and at bottom of the convective envelope for the 6.5 M 0 , Z = 0.012 model. The mass 
of the H-exhausted core shows evidence of efficient TDU whereas the He-shcll temperature 
shows that the majority of the TPs reach peak temperatures greater than 300 x 10 6 K, the 
temperature at which the 22 Ne neutron source is activated. HBB prevented several models 
in Table 1 from becoming carbon rich, and also resulted in low 12 C/ 13 C ratios that approach 
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the equilibrium ratio of ~ 4. The He/H and N/O ratios are also higher for models with 
HBB compared to the lower mass models. HBB is more efficient at lower metallicities for a 
given mass owing to the higher temperatures reached at the base of the convective envelope. 
For example, the 6.5 M 0 , Z = 0.012 model (corresponding to the Asplund et al. solar 
abundances) displays a higher N/O ratio than the Z = 0.02 model of the same mass. The 
final N abundance is similar in each of these two models, despite a difference of a factor 
of 1.86 in the initial N abundances, showing that CNO processing is more efficient in the 
low-metallicity case. 


3.2. The nuclear network 

The most important addition to our models for this study is the extension of the nu- 
clear network to include elements heavier than iron up to niobium (M. van Raai et ah, in 
preparation). We have also made an update of our reaction-rate library - starting from the 
library described in Karakas et ah (2007) we have included neutron capture cross sections 
from the Bao et ah (2000) compilation. We include 166 species from protons to sulfur and 
iron through to Nb, with a total of 1285 reaction rates corresponding to all /5-decay, p, a, 
and n-capture reactions on all species in the network. We also have two other networks 
that we use to obtain information for the S and Ar isotopes. The first has 125 species and 
includes all stable species from H to 62 Ni, while the second has 156 species, and includes all 
stable species from H to arsenic (Karakas et ah 2007). 

As done in Karakas et ah (2007), we use a “double neutron-sink” description (Herwig, 
Langer, & Lugaro 2003) to account for neutron captures on species heavier than that at the 
end of the network ( 9 'Nb in this case). The two artificial species are linked by the reactions 
97 Nb(n, y) 98 g and 98 g(n, L) 98 g, where 98 g has an initial abundance equal to the sum of the 
solar abundances from Mo to Bi. The second artificial particle, L , is equivalent to the number 
of neutrons captured beyond Nb. The ratio (L/ 98 g) is a description of the neutrons captured 
per seed nucleus and could in principle be related to the s-process distribution. Note that 
in the 156 and 166 networks, the species 34 S is the sum of all species from 34 S to Mn, and 
the reaction 34 S(n, 7 ) 35 S is assigned an averaged cross section value in order to represent all 
nuclei from 34 S to Mn. 
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3.3. Neutron sources in intermediate-mass AGB stars 

The s process is driven by the production of free neutrons that are subsequently captured 
by Fe-peak seed nuclei to form heavier elements. The 22 Ne(a, n) 25 Mg reaction was first 
suggested as a neutron source in stars by Cameron (1960). Later, Truran & Iben (1977) 
and Cosner et al. (1980) suggested it was the dominant source in intermediate-mass AGB 
stars, because the He-shclls of these stars reach high enough temperatures ( T > 300 x 
10 6 K) to allow for this reaction to be efficiently activated (see the middle panel of Fig. 3 
as an example). These temperatures are reached only in the last few TPs of a lower mass 
star. Truran & Iben (1977) found that the 22 Ne source results in enhanced levels of the 
n-rich element Rb owing to the high neutron densities (up to ~ 10 13 neutrons/cm 3 ), as 
well as increases in the heavy Mg isotopes 25 Mg and 26 Mg from the competing reactions 
22 Ne(a, n) 25 Mg and 22 Ne(a, y) 26 Mg (Kaeppeler et al. 1994; Karakas et al. 2006). Fenner et al. 
(2003) compared results from a chemical evolution model with observations of the neutron- 
rich Mg isotopes and concluded that an extra production site besides Type II supernovae 
was necessary, and that massive AGB stars are a good candidate for such a site. There is, 
however, some uncertainty at what Galactic epoch AGB stars started contributing to the 
chemical enrichment of the Galaxy (Simmerer et al. 2004; Melendez & Cohen 2007). It is 
also unknown to what extent intermediate-mass stars contribute to the Galactic inventory 
of s-process elements (e.g. Travaglio et al. 2004). 

The other potential source of neutrons in AGB stars is the 13 C(a,n) 16 0 reaction, which 
operates at lower temperatures (T > 90 x 10 6 ) than the 22 Ne source. Observational and 
theoretical evidence suggests this is the dominant neutron source in low-mass AGB stars 
(Smith et al. 1987; Gallino et al. 1998). To operate efficiently, this reaction requires more 
13 C than is left over from CN cycling in the H-shcll; hence some mechanism to mix protons 
from the H-rich envelope into the intershell is needed to produce the extra 13 C. In our 
models, protons are mixed into the intershell region by artificially adding a partial mixing 
zone (PMZ) at the deepest extent of each TDLI. These protons are quickly captured by 
the abundant 12 C to form 13 C and 14 N, resulting in the formation of a 13 C pocket. In the 
13 C pocket, neutrons are liberated by the reaction 13 C(a, n) 16 0 during the interpulse period 
(Straniero et al. 1995), in contrast to the 22 Ne source, which operates during TPs. The 
timescales for neutron production during the interpulse are much longer (> 10 3 years) than 
during the convective pulse (~ 10 years), resulting in much lower neutron densities (~ 10 7 
neutrons/cm 3 ) than the 22 Ne source. Together, the timescale for neutron production and 
the neutron source determine the resulting s-process element distribution. The details of 
how the 13 C pocket forms and its extent in mass in the He- intershell are still unknown, 
although various mechanisms have been proposed, including convective overshoot, rotation, 
and gravity waves; see Herwig (2005) for a discussion of the relative merits of each mechanism. 
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The importance of the 4 * * * * * * * * 13 C pocket and the 13 C(cpn) 16 0 reaction for nucleosynthesis 
in intermediate-mass stars is neither clear nor well studied. Adding a PMZ artificially into 
massive AGB models has little effect on the nucleosynthesis of nuclei lighter than Fe (Karakas 
et al. 2006) and on the production of Ge in 5 M 0 , Z = 0.02 models (Karakas et al. 2007). The 
first point can be understood by noting that the 13 C pocket is about ~ 10% of the mass of the 
He-intershcll. Since the He-intershcll is already smaller in massive AGB stars by about one 
order of magnitude compared to lower-mass stars, the inclusion of the PMZ has little effect 
on the surface abundances of light elements, which are relatively abundant in comparison to 
heavy elements produced via the s process. The effect of the 13 C pocket on the production 
of s nuclei is less clear, given their low initial abundances. Se and Kr elemental abundances 
derived from Type 1 PN spectra, in combination with observations of the elements Rb and 
Zr in massive AGB stars in the Galaxy, LMC, and SMC, may provide a way of distinguishing 
between the relative importance of the two neutron sources in intermediate-mass AGB stars. 

For the 3, 4, 5, and 6.5 M 0 solar composition models we present results with and without 
a PMZ (see Table 2 for the PMZ masses used in each model). No PMZ was included in any 
of the low-metallicity models, with the exception of the 2.5 M 0 , Z = 0.008 model, where we 
only include results with a PMZ. The PMZ masses were chosen such that the resulting 13 C 
pockets were ~ 10% to 15% of the mass of the He- intershell, consistent with the 13 C pocket 
masses used by Gallino et al. (1998). The method we use to include a PMZ is described 
in Lugaro et al. (2004), and is similar to that used by Goriely & Mowlavi (2000). At the 
deepest extent of each TDU episode, we add an exponentially decaying proton profile that 
covers the mass specified in Table 2. This method is different from that used by Gallino 
et al. (1998), who include a 13 C profile directly into the intershell. 


4. Model Results 

4.1. Germanium to Strontium 

Surface abundances at the tip of the TP-AGB for Zn to Sr are presented in Table 2, 

and Figure 4 depicts enrichments of Se and Kr as a function of pulse number in selected 

models. At solar metallicity, the 3M 0 models with a 13 C pocket produce the largest surface 

enrichments of s-process nuclei, with final [Se/Fe] and [Sr/Fe] abundances of ~ 0.50 and 0.85 

for the Z = 0.02 model, and 0.56 and 1.2 for the Z = 0.012 model. In the lower metallicity 

models, the largest surface enrichments are seen in the 2.5 M 0 , Z = 0.008 model with final 

[Se/Fe] and [Sr/Fe] abundances of 0.60 and 1.3, and in the 5M 0 , Z = 0.004 model with 

0.90 and 0.93, respectively. 


r^j 
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The solar metallicity intermediate-mass (M > 5M 0 ) models exhibit little s-process 
enrichment ([Kr/Fe] and [Sr/Fe] < 0.25). Using Asplund et al. (2005) solar abundances, 
the 6.5 M 0 , Z = 0.012 model shows slightly larger enhancements ([Kr/Fe] = 0.40, [Sr/Fe] 
= 0.30), albeit only when a PMZ is added. The inclusion of a 13 C pocket has little effect 
(< 0.1 dex) on the production of Se and Ge, as noted by Karakas et al. (2007), but it does 
enhance the Kr and Sr abundances. Only the lower-metallicity 5 and 6M 0 models produce 
substantial enrichments of s-process elements. These models experience a large number of 
TPs (> 60), each with efficient TDU, resulting in significant surface enrichments. No 13 C 
pocket was used for these models, and thus it is clear that the higher temperatures obtained 
during TPs efficiently activates the 22 Ne neutron source. Moreover, the number of neutrons 
per seed nucleus is larger at lower metallicity owing to the smaller abundance of seed nuclei 
(e.g., 56 Fe). 

We find that Kr is enhanced in the stellar envelopes to a greater degree than Se, in 
agreement with the observational results of Sterling et al. (2007) and SD08. This is because 
Kr lies nearer to the first s-process peak than Se, and has an isotope ( 86 Kr) with a magic 
number of neutrons 3 . Furthermore, the larger s-process enhancement of Kr relative to Se 
in the AGB models are in agreement with the chemical evolution results of Travaglio et al. 
(2004), who find that the main s-process contributions of AGB stars to the Solar System Se 
and Kr abundances are 14% and 30%, respectively (the weak s-process, operative in massive 
stars, contributes 25% and 20% to the solar Se and Kr abundances, while the r-process 
provides the remainder). 

The [Ge/Fe] abundances in Table 2 are different from those given in Table 2 in Karakas 
et al. (2007), owing to small changes to the neutron-capture cross sections used in the 166 
and 156 species networks. There is a 0.08 dex change in the [Ge/Fe] abundance for the 3 M 0 , 
Z = 0.02 model with PMZ of O.OO2M 0 . This difference of 0.08 dex is less than found when 
considering other model uncertainties. For example, there is a 0.16 dex difference in the 
[Ge/Fe] abundances from the two different 5M 0 , Z = 0.02 models with different mass-loss 
laws. 

Efficient activation of the 22 Ne neutron source leads to positive <5 86 Kr/ 82 Kr ratios 4 , as 
can be seen in Table 2 for the most massive AGB models. These positive <5 86 Kr/ 82 Kr ratios 


3 Nuclei with a magic number of neutrons are relatively stable against n capture due to their small n- 
capture cross sections, and are the cause of the peaks in the s-process enrichment distribution. 

4 The 5 notation is commonly used in meteoritics (see, e.g., Lugaro et al. 2003a), and is a variation from 
solar in parts permil, or parts per thousand (where a 100 permil change in the isotopic value is equivalent 
to a 10% variation). 
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are a result of neutron densities in excess of 10 8 cm -3 , which enables neutron captures on 85 Kr 
that produce 86 Kr. If the density is lower, as is the case when 13 C(a,n) 16 0 is the neutron 
source, the branching is bypassed and 85 Rb is produced instead, resulting in negative values 
of the <5 86 Kr/ 82 Kr ratio. The Z = 0.02 models reflect this trend, showing an increase in the 
5 86 Kr/ 82 Kr ratio with mass, with 5 M 0 the minimum mass for efficient activation of the 22 Ne 
source. The ratio is negative in the 3-4 M 0 models, where the 13 C reaction dominates as the 
neutron source and the 22 Ne reaction is only marginally activated. Values of 5( 86 Kr/ 82 Kr) 
measured in meteoritic stellar SiC grains from AGB stars range roughly from —500 to 900, 
and increase with the grain size (Lewis et al. 1994). Negative values are well explained by 
AGB model predictions, however, it is difficult to explain the large observed positive values 
because they are only produced by HBB models, where the C>0 condition necessary for the 
formation of SiC grains is not attained, except in the lowest metallicity, Z = 0.004 models. 
Further detailed work is required on this issue. 

The resulting elemental abundance pattern also differs, as Sr is produced in favor of Rb 
in the low-neutron density case, resulting in low Rb/Sr ratios (~ 0.05). This is consistent 
with observations of low- mass AGB stars (Lambert et al. 1995; Abia et al. 2001). In the 
high-neutron density case, Rb is produced instead, as observed in massive O-rich AGB stars 
(Garcia-Hernandez et al. 2006). Preliminary results for Zr and Rb from intermediate-mass 
AGB stars are discussed in van Raai et al. (2008), and will be presented in more detail in a 
forth-coming paper (M. A. van Raai et ah, in preparation). 


4.2. The iron-peak elements Iron and Zinc 

In Table 2, we show the [Zn/Fe] ratios and AFe values for each model, where AFe is the 
logarithmic change to the surface iron abundance from the main sequence to the tip of the 
TP- AGB. The AFe values are all very nearly zero, indicating that the surface abundance 
of this element is unaltered by AGB nucleosynthesis. The maximum effect among the Z = 
0.02 models is 2%, and even for the Z = 0.004 case the Fe abundance decreases by less than 
5%. This makes Fe a suitable reference element for estimating metallicities in AGB stars. 

Although the surface Fe abundance is not significantly changed by nucleosynthesis dur- 
ing the AGB, the abundance of Fe in the He-intershell can be depleted by up to ~ 1 dex in a 
massive AGB star (e.g., 6.5 M 0 model). In the intershell 56 Fe is destroyed by nentron cap- 
tures, and as the dominant isotope it primarily determines the elemental abundance (even 
if the abundances of the rare neutron- rich isotopes 5 'Fe and 58 Fe increase). In lower mass 
models (e.g., 3 M 0 ), depletion of Fe is restricted to the 13 C pocket that encompasses 10 to 
15% of the mass of the intershell. The next TP homogenizes this small, Fe-depleted region 
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with the much larger Fe abundance in the rest of the shell. Onr 3 M 0 results are consistent 
with those of Werner & Herwig (2006), who found that AGB models cannot account for 
the large Fe depictions found in some central stars of PNe, particularly the hot, H-deficient 
stars of the PG1159 class (Miksa et ah 2002; Rauch et al. 2008). Werner & Herwig (2006) 
estimated that these central stars have mean masses of ~ O.62M 0 , and thus represent the 
evolutionary remnants of lower-mass progenitors than those addressed here, and should be 
compared with models of stars with < 3M 0 . Note that similar or even more extreme Fe 
deficiencies - along with dramatic overabundances of n-capture elements - are seen in older 
white dwarfs where mechanisms such as gravitational settling and levitation by selective 
radiation pressure play important roles in establishing the surface abundances (Chayer et al. 
2005; Werner & Herwig 2006). 

In PNe, Fe can be depicted into dust by large factors that may vary from object to 
object (Perinotto et al. 1999), and even in different regions within the same nebula (Sterling 
et al. 2005), so that its gas-phase abundance is not a reliable indicator of the true elemental 
abundance. The heavier iron-group element Zn exhibits much milder depletions in the ISM, 
particularly in the warm phases (e.g., Welty et al. 1999). Consequently, Zn is frequently 
used as a proxy for iron and the iron-group, for example in damped Lyman-a absorbers 
(e.g., Akerman et al. 2005). However, the use of Zn as an indicator of the initial composition 
of an AGB star or PN requires that the change in the abundance of Zn due to internal 
nucleosynthesis is small compared to its initial abundance in the star. It can be seen from 
Table 2 that the predicted changes in the Zn abundance are modest or negligible for most of 
the models. This is consistent with the conclusions of Travaglio et al. (2004), who estimated 
the total contribution of AGB stars to the solar Zn to be ~ 3%, and with Cayrcl et al. 
(2004) ’s finding that Zn production in the Milky Way is dominated by other sources (i.e., 
explosive Si burning in massive stars). However, at low metallicities, the models show a large 
enough increase that the Zn enrichment may be detectable. 

These results can be compared to the Zn abundances derived for two Type I PNe 
(DGS08). Of the models presented here, the Galactic disk PN M 1-40 is most closely matched 
by the 6.0 M 0 model with Z = 0.02, given its solar Ne and S abundances and N/O = 1.13 
(Gorny et al. 2004). This model predicts a very small enrichment of Zn, [Zn/Fe] ~ 0.038, 
which is consistent with the measured value (see § 2), even allowing for a possible small 
depletion effect. The case of LMC SMP 62 is more problematical. According to its a-elcment 
abundances, the most appropriate models are those for Z = 0.008. At this metallicity, all 
three of the models in Table 2, calculated for 2.5, 5.0, and 6.0 M 0 respectively, predict 
similar Zn enhancements of ~0. 15-0.2 dex. However, as described in § 2, the measured Zn 
abundance in this PN is lower, rather than higher, than those of the a-elements (e.g., Ne, S), 
in contradiction to the models. The disagreement is even worse for Se, which is predicted to 
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be enhanced by 0.4-0. 6 dex, depending on the initial mass. A modest enrichment in both Zn 
and Se is possible, if their initial abundances were low enough relative to Fe, but nevertheless 
it remains difficult to satisfy the observational constraints for LMC SMP 62 with the listed 
models. We note, however, that these particular models also predict N/O values either lower 
(the 2.5 M 0 model) or higher (the 5 and 6 M 0 models) than the observed value for this 
PN; in § 5 we discuss other possible scenarios for explaining the observed properties of LMC 
SMP 62. 


4.3. Neon through to Argon 

Abundances of S, Cl, and Ar in PNe are powerful tools for studying the chemical evolu- 
tion of galaxies and stellar systems, under the assumptions that these elements are unaltered 
by AGB evolution (Dopita et al. 1997; Stasihska et al. 1998; Leisy & Dennefeld 2006; Idiart 
et al. 2007) and are produced in the same relative proportions in massive stars (Stasihska 
et al. 1998). The case of Ne is not as clear, although the remarkably constant Ne/O ratio 
of ~ 0.20 (by mass) observed in PNe with a wide variety of chemical compositions indicates 
that in most cases this ratio is not altered during the AGB. 20 Ne, the most abundant isotope 
of Ne, is not altered significantly by H or He-burning in AGB stars, whereas 22 Ne can be 
produced by successive a-captures onto 14 N. Therefore an enhancement in the Ne/O ratio 
is expected only when 22 Ne is sufficiently enriched that its abundance is comparable to or 
exceeds that of 20 Ne. Karakas & Lattanzio (2003b) predicted that such an increase would 
only occur in a narrow range of progenitor masses, ~ 2.5 to ~ 3.5 M 0 depending on rnetal- 
licity. Our current models are generally consistent with these earlier results, although from 
Table 3 it is clear that the 6.5 M 0 , Z = 0.012 model and the 5 and 6 M & , Z = 0.008 and 
0.004 models have enhanced Ne/O ratios as a result of some O destruction by HBB, and 
mild enrichments of Ne via dredge-up (maximum Ne increase of ~ 0.3 dex for the 5, 6M 0 
metal-poor models). Therefore, Type I PNe are not expected to be significantly enriched in 
Ne except at low metallicities. 

The 22 Ne is produced primarily in the TPs themselves, but in lower mass models (2.5 and 
3 M 0 ), the partial mixing of protons can further enhance the Ne/O ratio. This is because the 
partially mixed zone results in a pocket of 13 C and 14 N, thus increasing the 14 N abundance in 
the intershell to values higher than that from CN-cycling during the interpulse. The next TP 
quickly converts the 14 N into 22 Ne. The resulting 22 Ne abundances after the TP are slightly 
higher than in the case with no partially-mixed zone. After many TPs and TDLI, the PMZ 
results in an increase in the elemental Ne abundance. The observations by Bernard-Salas 
et al. (2008) are intriguing in this respect. They determined Ne/S ratios in a sample of 25 
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PNe in the Magellanic Clouds, and found that four of the objects show high Ne/S ratios. 
For two of the objects, they attribute the high ratios to low S abundances, but in the other 
two PNe, the high Ne/S values could be due to Ne enrichments. 

Table 3 shows the surface abundance predictions for S, Cl, and Ar. These elements 
generally are not altered by AGB nucleosynthesis processes by more than the observational 
errors (~ ±0.1 — 0.15 dex), consistent with the observations of PNe by Stasinska et al. (1998). 
Neutron captures in the He-shell can result in variations in the isotopic ratios involving the 
stable Cl and Ar isotopes; this is demonstrated in the last column of Table 3 where we show 
the permil variations for 38 Ar/ 36 Ar. Values of <5( 38 Ar/ 36 Ar) measured in stellar SiC grains 
range from zero (solar composition) to positive values up to ~ 200 (Lewis et al. 1994), in 
qualitative agreement with our model predictions. Kahane et al. (2000) compared AGB 
model predictions to measurements of Cl, Mg, Si, and S isotopic ratios in the cool carbon 
star IRC±10216, and inferred an initial mass of < 2 M 0 based on the accuracy of the Cl and 
Mg data. Kahane et al. concluded that a larger mass progenitor would lead to too high Cl 
and Mg isotopic ratios at the surface, although our models suggest that progenitor masses 
as high as ~ 3 M 0 produce a good fit to the observed data, within the errors quoted for the 
isotopic ratios. In the case of our low-metallicity massive AGB models, neutron captures 
result in modest increases (~ 0.25 dex) in the elemental Cl surface abundance. 

The phosphorous abundance can be determined in hot post- AGB objects such as the PG 
1159 stars, whose envelopes are composed of He- intershell material. Observations indicate 
that PG 1159 stars have solar P abundances, whereas 3 M 0 models predict P abundances 
ranging from four times solar to ~ 25 times solar in the He-intershell, depending on the 
amount of extra mixing included in the calculations to produce the 13 C pocket (Werner 
& Herwig 2006). Our He-shell predictions are consistent with Herwig’s model, in that we 
estimate P enhancements from 1.3 to ~ 9 times solar for the 3 M 0 model, without and with 
a PMZ respectively. The large P abundances in the He-intershell abundances do not result in 
large surface enrichments (< 0.1 dex in all solar metallicity models), so that if measurements 
could be made of the P abundances in PNe, we would not expect to see variations from solar. 

Sulfur is also observed in PG 1159 stars, with abundances varying from 0.01 times solar 
to solar (Werner & Herwig 2006). It is not clear if this range reflects the initial abundance 
of the star, various degrees of nuclear processing, condensation of S into dust or problems 
in determining accurate S abundances. Note that S abundances in several Galactic (Henry 
et al. 2004) and some Magellanic Cloud PNe (Bernard-Salas et al. 2008) are low relative to 
O. While infrared observations of S 3+ in Galactic PNe seem to indicate that uncertainties in 
corrections for unobserved ionization stages are not the reason for the low S abundances, it is 
possible that S may be depleted into dust in some PNe (especially C-rich objects; Pottasch & 
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Bernard-Salas 2006). AGB models do not predict large depletions of S in the H e-shell, with 
S abundances varying between 0.6 to 0.9 times solar (Werner & Herwig 2006). Our values 
lie from no change to ~ 0.9 times solar. The results for P and S are puzzling, indicating on 
the one hand that a smaller PMZ is required (to produce less P), or, on the other hand that 
more mixing is required (to destroy more S). 


4.4. Oxygen in intermediate-mass AGB stars 

Oxygen constitutes a large fraction of the total metal content of a star, and in this sense is 
a more fundamental measure of metallicity than Fe/H, which historically has been used as the 
primary indicator of metallicity due to its ease of measurement. While it can be challenging to 
determine O abundances in stars (e.g., Ramirez et al. 2007), O/H is easily measured in ionized 
nebulae because of its bright emission lines, and the resulting abundances are considered 
to be among the most accurate for any element. According to the current paradigm for 
galactic chemical evolution, O production is dominated by massive stars (M > 10 M 0 ). 
Nucleosynthesis models of AGB stars indicate that, collectively, these stars do not contribute 
significantly to the production of O (Forestini & Charbonnel 1997; Karakas & Lattanzio 
2007), although there is evidence for dredge-up of self-produced oxygen in a few special 
cases, including PNe from metal-poor progenitor stars (e.g., Pequignot et al. 2000; Dinerstein 
et al. 2003). The inclusion of diffusive convective overshoot at the inner edge of the flash- 
driven convective zone during a TP can lead to O abundances in the intershell of ~ 20% 
(Herwig 2000), in contrast to standard models that predict 2% by mass or less (Boothroyd & 
Sackmann 1988; Karakas et al. 2002). This can lead to increases in the surface composition 
of O, particularly at low metallicities. Evidence for this overshoot comes from PG 1159 
stars that show compositions consistent with Herwig’s models (Werner & Herwig 2006). 
It is unclear however, if this overshooting occurs for all AGB stars or only for AGB stars 
that produce PG 1159 stars as a consequence of late and very late TPs. For example, the 
observed abundance analyses in intrinsic and extrinsic AGB stars, as well as in SiC grains, 
suggest that such an overshoot is uncommon (Lugaro et al. 2003b). This is because in the 
region of the He-rich intershell where primary O is synthesized, the temperature is large 
enough for efficient activation of the 22 Ne(a, n) 25 Mg neutron source, even in low-mass AGB 
stars. In this case, the resulting neutron density is very large (up to ~ 10 12 neutrons/cm 3 ) 
in contrast to observations of low-mass AGB stars that indicates neutron densities of about 
~ 10 7 neutrons/cm 3 (Abia et al. 2001). 

Intermediate-mass AGB stars on the other hand, destroy O via HBB when the tem- 
perature at the base of the convective envelope is sufficiently high (T > 80 x 10 6 K). The 



O-Na anti-correlation observed in globular clusters stars has in fact been associated with 
HBB in intermediate-mass AGB stars, because HBB can qualitatively account for the de- 
struction of O along with the production of Na. Moreover, SD08 noted that the Ar/O, S/O, 
and Cl/O ratios of Type 1 PNe are approximately a factor of two larger on average than in 
non-Type 1 objects. Since Ar, S, and Cl are not significantly affected by nucleosynthesis in 
low- or intermediate-mass stars (see § 4.3), they interpreted this as evidence for O destruc- 
tion during HBB in Type 1 PN progenitors. This result had also been found in previous 
investigations (Marigo et al. 2003; Leisy & Dennefeld 2006; Pottasch & Bernard-Salas 2006). 

The amount of O destroyed by HBB depends on the temperature at the base of the 
convective envelope. The temperature, in turn, depends on a number of factors including the 
initial mass and metallicity, with hotter temperatures found in models of increasing mass for 
a given Z, or in models of decreasing metallicity, at a given mass. In addition, the predicted 
efficiency of HBB is sensitive to the treatment of convection and mass loss in the models, 
neither of which are well- constrained empirically. In particular, the convective efficiency 
affects the luminosity and the mass- loss rate (Ventura & D’Antona 2005), and hence these 
properties are sensitive to the convective model that is utilized. Mass loss determines the 
length of the HBB phase, since HBB does not operate once the envelope mass drops below 
about 1.5 M 0 . 

In Table 3 we present the O abundances at the surface of the model star at the tip 
of the AGB. The O abundance remains unaltered in the lower mass models (3 and 4M 0 , 
Z = 0.02), whereas O destruction takes place in models with HBB. Out of the Z = 0.02 
models, the largest amount of O destruction is seen at 6.5 M 0 which has a peak temperature 
at the base of the envelope of 86xl0 6 K. The effect of metallicity is seen most clearly by 
comparing the 6.5 M & models with Z = 0.02 and Z = 0.012, with the latter model computed 
assuming Asplund et al. (2005) C, N, O initial abundances. In the Z = 0.012 model the peak 
temperature is higher, resulting in 0.15 dex destruction of O compared to only 0.09 dex in 
the Z = 0.02 model. Note that the O that is destroyed is all converted into N. The largest 
O destruction of 0.26 dex occurs in the 6M 0 , Z = 0.008 model. Overall, none of the models 
suffered O destruction of more than 0.3 dex. 

It is not clear that the low O abundances relative to S, Ar, and Cl in Type I PNe can be 
attributed to HBB. While modeling uncertainties in the treatment of convection and mass- 
loss could in principle allow for more efficient O destruction during HBB than we report, this 
inevitably leads to large increases in the N/O ratio. This result agrees with the findings of 
Marigo et al. (2003), whose synthetic models with efficient ON-cycling produced N/O ratios 
considerably higher than those observed in Type I PNe. Those authors suggested that the 
low O abundances may be explained if the initial metallicities of Type I PN progenitors are 
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low. However, this does not explain the approximately solar abundances of S and Ar in these 
objects. Moreover, our models constrain the metallicities of Type 1 PN progenitors to be 
Z > 0.4 Z & in order to reproduce the observed Se and Kr abundances. This important issue 
remains unresolved. 


5. Discussion 

In Fig. 5 we illustrate the s-process abundance predictions from a selection of the stellar 
models. The gray-hatched box represents the approximate region occupied by Type I PNe 
s-process abundances. Typical observational uncertainties allow the maximum enrichment 
to reach « 0.6 dex in some cases. The intermediate-mass AGB models are, in general, a 
good match to the observed abundances of Se and Kr in Type I PNe. Models in Table 2 
that show enhancements of Se and Kr greater than 0.5 dex are the 2.5 M 0 model, the 3 M & 
models, the 6 M 0 , Z = 0.008 model, and the low-metallicity 5 M 0 , Z = 0.004 model, which 
gives large enhancements for all elements between Zn and Sr. The results for 3M 0 with 
a PMZ are only in borderline agreement with the observations, although the results with 
no PMZ are consistent with no s-process enrichments. In either case the final composition 
is carbon-rich, indicating that this model is unlikely to be the precursor of a Type I PN. 
Indeed, the results for the 3M 0 models are closer to that found for the 2.5 M & , Z = 0.008 
model which has a distinct, non-Type I PNe composition, with a high final C/O ratio, and 
considerable overabundances of all s-process elements (> 0.5 dex for Se to Sr). The results 
for the 5M 0 , Z = 0.004 model suggests that Type I PNe in the Galaxy do not descend 
from low-metallicity stars, as suggested, for example, by Marigo et al. (2003), and that a 
minimum metallicity of ~ 0.4 Z Q is required for these PNe. This result is not unexpected, 
given that Type I PNe evolve from more massive stars than other PNe, and hence trace a 
younger population. On the other hand, these model results suggest that Type I PNe in the 
Magellanic Clouds could be more highly enriched in light n-capture elements such as Se and 
Kr than their Galactic counterparts. This is a strong motivation for further observations 
of s-process elements in PNe belonging to the Magellanic Clouds or to other metal-poor 
populations. 

Observations of other trans-iron elements in PNe in addition to Se and Kr could be 
used to constrain the mass and metallicity ranges of Type I PNe. For example, Sr and 
Rb abundances are potentially useful discriminants of progenitor mass. Sr is considerably 
enriched in the 3M 0 models ([Sr/Fe] ~ 0.85 — 1.23), but is slightly if at all enhanced in the 
6 and 6.5 M 0 models ([Sr/Fe] < 0.4 dex). Note that the inclusion of a PMZ in the 6.5 M 0 , 
Z = 0.012 model increases the level of Sr enrichment expected, from 0.05 to 0.3 dex. 



Of these next two elements heavier than Kr (element number 36), Rb (element 37) is 
more likely to be measurable in a meaningful way in nebulae than Sr (element 38), for the 
following reasons. First, Sr is refractory and can be strongly depicted into dust grains. In 
fact, based on its condensation temperature (Lodders 2003), Sr is expected to be incorporated 
into dust to a greater degree than Fe, whose gaseous abundance is depleted by 1-2 orders 
of magnitude in PNe (Perinotto et al. 1999; Sterling et al. 2005) and the ISM (Welty et ah 
1999). Therefore, measurements of the gaseous Sr abundance in PNe are not likely to reveal 
useful information about s-process enrichments. Second, the only emission lines of Sr in the 
optical spectral region originate from relatively highly ionized states that are not expected to 
be abundant, except in a few of the highest-excitation PNe. Sr emission has been identified 
in only one PN to date (Pequignot & Baluteau 1994; Zhang & Liu 2005). On the other 
hand, based on its condensation temperature (Lodders 2003) Rb does not deplete into dust 
as readily as Sr, and thus its gaseous abundance should be more representative of the overall 
elemental abundance. Rb emission has been identified in a number of PNe (Pequignot & 
Baluteau 1994; Zhang & Liu 2005; Sharpee et ah 2007) and should be detectable in PNe with 
a broader range of ionization than Sr. This discussion suggests that Rb is a useful element 
with which to constrain stellar models to PN observations. 

Furthermore, Rb is observed to be substantially enhanced in massive Galactic OH/IR 
stars with the maximum enhancement [Rb/Fe] ~ 2 ± 1 dex (Garcia- Hernandez et ah 2006). 
van Raai et ah (2008) have made a preliminary comparison between Zr, Rb, and Li predic- 
tions from 5 to 6.5 M 0 solar composition models to the observations of the OH/IR stars. The 
best match between stellar models and the observations was obtained with the 6.5 M 0 model 
computed by 1) not including a 13 C pocket, and 2) by accounting for extra TPs not modeled 
in detail, as done by Karakas et ah (2007). The extra TPs and TDU mixing episodes may 
occur because not all the envelope mass was lost when convergence difficulties ended the 
computation. The main uncertainty is the behavior of the TDU efficiency at small envelope 
masses, with studies generally finding a decrease in the efficiency with decreasing envelope 
mass (Straniero et ah 1997; Karakas et ah 2002). Further studies on the evolution of the 
TDU efficiency with decreasing envelope mass are needed to help settle this issue. The tab- 
ulated results presented in this study were calculated without the inclusion of these extra 
TPs and can be considered lower limits to the final PN abundances, noting that it is highly 
uncertain how much of an impact the final TPs would have on the the final PN abundances. 

The final few TPs can potentially have a significant effect on the final abundance of the 
star and its resulting PN. For example, Frost et ah (1998) noted that intermediate-mass AGB 
stars may become luminous, optically obscured carbon stars near the end of the TP-AGB, 
when mass loss has removed much of the envelope, extinguishing HBB but allowing dredge- 
up to continue. Under these conditions the envelope mass is signficantly reduced, resulting 
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in only minimal dilution of the intershell material that is dredged up. This can allow for 
large enhancements to the surface composition of C and other He-burning products. In the 
models with HBB included in this study, all experienced a few TPs and TDU episodes after 
the cessation of HBB, with some of the lower metallicity models becoming carbon rich (see 
Table 1, i.e., 6 M e , Z = 0.008; 5 and 6 M 0 , Z = 0.004) before convergence difficulties ended 
the computation. Note that it is after this stage that we consider the effect of remaining 
TPs not modeled in detail using a synthetic AGB algorithm. 

In van Raai et al. (2008) the inclusion of these remaining TPs only had a large impact on 
the surface abundances of the 6.5 M 0 , Z = 0.02 model (out of the 5, 6 and 6.5 M 0 , Z = 0.02 
models), and only when assuming that efficient TDU occurs at a very small envelope mass 5 . 
A final [Rb/Fe] ~ 0.9 was obtained in the 6.5 M 0 , Z = 0.02 model, in agreement with the 
observed Rb abundances of OH/IR stars given the large uncertainties. Along with increases 
in Rb, the [Se, Kr/Fe] abundances were increased to 0.6 and 0.7 dex (respectively), above 
the observed Type I PN limit of 0.3 dex. This model was estimated to have 7 remaining 
TPs, with the last TP occurring with a very small envelope mass of ~ O.O3M 0 . The small 
envelope mass during the last TDU leads to the situation discussed above: A substantial 
increase in the surface abundances, caused by the mass of the envelope being only ~ 10 
times greater than the mass of the intershcll. It is questionable whether such an efficient 
TDU would occur for such a low envelope mass; so if we exclude this last mixing episode, the 
[Se, Kr/Fe] abundances are ~ 0.4 and 0.45, respectively, in agreement with observed Type I 
PN abundances. This analysis may indicate that the OH/IR stars are not the progenitors 
of most Type I PNe and that there are possible observational biases in the sampling of each 
population. That is, the OH/IR population may sample the most massive AGB stars (6- 
7 M 0 ), with the largest [Rb/Fe] (and [Kr/Fe]) abundances, while the Type I PN population 
may come from slightly less massive AGB stars, with lower Kr and Rb abundances. There is 
also the problem of explaining the large Rb abundances in the OH /IR stars, given that they 
seem to require the TDU to operate efficiently at small envelope masses. More observations 
of Rb in OH/IR stars and PNe will help to settle these issues. 

Modeling uncertainties can dramatically affect AGB nucleosynthesis predictions. As 
noted previously, different convective models have a large impact on the nucleosynthesis of 
light elements (e.g., N, O, Na) in AGB models with HBB. Improvements in our understanding 
of convection are desperately needed but this may only come about through multidimensional 
studies (e.g., Dearborn et al. 2006; Herwig et al. 2006; Meakin & Arnett 2007). The mass-loss 
law used in the computations is similarly important, because this can determine (for a fixed 


5 Besides the inclusion of the extra TP the models presented here and in van Raai et al. (2008) are the 


same. 
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convective model) the number of TPs and mixing episodes. For example, if we artificially cut 
the number of TPs by half for the 5 M 0 , Z = 0.004 model (40 TPs instead of the computed 
80), the final [Kr/Fe] ratio is ~ 0.35 dex instead of 1.17 dex. This simple analysis indicates 
that we need better constraints on the mass loss from low-metallicity, massive AGB stars if 
we are to provide constraints on the metallicity range of Type 1 PNe. 

Another important uncertainty relates to the unknown formation mechanism for 13 C 
pockets in AGB stars (Gallino et al. 1998; Herwig 2005). These pockets are necessary to 
provide neutrons in low-mass stars. It is unclear, however, what role the 13 C neutron source 
plays in more massive AGB stars. We had hoped to ascertain the dominant neutron source 
from comparisons of Type I PN abundances to massive AGB stellar models, but the Se and 
Kr abundances alone are not sufficient. Selenium abundances were not strongly affected by 
the inclusion of a 13 C pocket in the 5 and 6.5 M 0 models, whereas the abundance of Kr was 
affected but not increased above 0.5 dex. Note that in three Type I PNe of SD08 where 
both Se and Kr were detected, an abundance of [Kr/Se] = 0.5 was found, providing some 
support for the inclusion of a 13 C pocket in massive AGB models; this conclusion is, however, 
limited by small number statistics. We predict that PNe should have larger enrichments of 
Kr than Se. While this result is valid regardless of whether or not a 13 C pocket is included, 
larger [Kr/Se] values are predicted in models with a PMZ (see Table 2). Firmer conclusions 
regarding the neutron source operating in intermediate-mass AGB stars may come from large 
data sets of many n-capture elements, including Br, Rb, and Xe. 

Finally, models with extra mixing may be necessary to reproduce the observed abun- 
dances of some PNe. For example, the composition determined for LMC SMP 62 (see §2 
and §4.2) suggests that it is a Type I PN, but its low Se and Zn abundances are at odds 
with our predicted abundances for intermediate-mass models. A 3 M 0 , Z = 0.02 model with 
rotation by Cantiello & Langer (private communication) had a N/O ratio of 0.43 after the 
first dredge up, compared to 0.32 from our model with no rotation, whereas the N/O value 
of LMC SMP 62 is ~ 0.4 - 0.5. During core H-burning, rotation mixes processed material to 
layers further out than they would otherwise be found. Rapid rotation may also inhibit the 
production of s-process elements during the AGB phase (Herwig et al. 2003), consistent with 
the lack of Zn and Se enrichments found in this object. An alternative is that LMC SMP 
62 evolved from a low- mass progenitor of less than 2 M 0 , which did not experience efficient 
TDU. Note that the low C/O ratio of 0.17 found for this object (Aller et al. 1987) is also 
consistent with a star that did not experience any TDLT mixing (suggesting an initial mass 
of < 1.25M 0 ). It is unclear if rotation can explain the high N/O ratio in such a low-mass 
star. Other extra-mixing processes (e.g., thermohaline mixing) may need to be invoked in 
this case. 



6. Conclusions 


Type I PNe have high He/H and N/O ratios, indicating that they are the descendants of 
intermediate-mass AGB stars with initial masses between ~ 3 to 8 M 0 . Sterling & Dinerstein 
(2008) found that Type I PNe exhibit significantly smaller enrichments of Se and Kr (< 
0.3 dex) on average than other PNe. We calculated s-process enrichments in a set of AGB 
models covering a range in mass from 2.5 to 6.5M 0 , and metallicity from 0.2 Z@ to Z & . The 
2.5 Mq, Z = 0.008 model was included to show an example of a model that would produce a 
clear non-Type I PN abundance signature, with [Se, Kr/Fe] > 0.5 and C/O ~ 4. The main 
conclusion is that the results for the 3-6. 5M 0 are a good match to the observed abundances. 
The only real exception out of the HBB models is the low-metallicity 5 M 0 , Z = 0.004 model 
that produced much larger enhancements in Se and Kr than observed. This suggests that 
Galactic Type I PNe do not descend from such low-metallicity objects. 

We also compare calculated abundances for selected intermediate-mass (Ne, P, S, Ar) 
and iron-peak (Fe, Zn) elements to observations of post-AGB stars and PNe. We find that 
the elemental abundances of P, S, Cl, Ar, Fe, and Zn are essentially unchanged by AGB 
nucleosynthesis, although there are isotopic shifts caused by neutron captures in the He- 
shell. These results justify use of elements such as S, Cl, and Ar in PNe as tracers of 
Galactic chemical evolution. 

It is difficult to reach firm conclusions about the neutron source operating in massive 
AGB stars from Se and Kr abundances in Type I PNe. Certainly it seems that at the least the 
22 Ne source, with efficient TDU, is required to produce enhancements in s-process elements 
in massive Type I PN progenitors. It is less clear whether a 13 C pocket is also required 
in these stars. Increases in the [Kr/Fe] abundance ratio were observed in the models when 
a 13 C pocket was included, but not beyond the amounts observed in Type I PN spectra. 
Obtaining abundances for more n-capture elements, particularly Rb and elements beyond 
the first s-process peak, from a large data set of OH/IR stars and PNe will help distinguish 
among the possibilities. 

Finally, only the models with HBB (M > 5 M 0 , depending on Z) show the high He/H 
and N/O ratios that define the Type I PN class. Given that 3M 0 stars are more common 
than 6.5M 0 stars (according to initial mass function and evolutionary time considerations), 
it may be necessary that another mixing process other than HBB is active in stars of the 
lower-mass range (3-4M 0 ), if these stars do in fact evolve into Type I PNe. Rapid stellar 
rotation in single ~ 3 M 0 stars may, for example, be able to account for the increased He and 
N abundances in the progenitor stars. The problems of convection and mixing, and of other 
modeling uncertainties (e.g., n-capture cross sections) requires further study in the context 
of s-process nucleosynthesis in intermediate-mass AGB stars. 
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Table 1: Details of the stellar models. 


Mass 

Z 

TPs 

T^max 

i He 

Tmax 
1 bee 

Mdred 

M env 

HBB? 

C/O 

12 C/ 13 C 

N/O 

He/H 

3.0 

0.02 

26 

302 

6.75 

8-l(- 

-2) 

0.676 

No 

1.27 

108 

0.303 

0.122 

4.0 

0.02 

17 

332 

22.7 

5.6(- 

-2) 

0.958 

No 

0.99 

76.5 

0.336 

0.118 

5.0 

0.02 

24 

352 

64.5 

5.0(- 

-2) 

1.500 

Yes 

0.77 

7.84 

0.542 

0.133 

5.0 

0.02 a 

37 

368 

57.0 

i-o(- 

-1) 

1.922 

Yes 

1.71 

15.5 

0.465 

0.139 

6.0 

0.02 

37 

369 

83.1 

5.8(- 

-2) 

1.791 

Yes 

0.38 

10.8 

1.251 

0.150 

6.5 

0.02 

40 

368 

86.5 

4.7(- 

-2) 

1.507 

Yes 

0.40 

11.6 

1.205 

0.154 

3.0 

0.012 

21 

307 

7.23 

9.2(- 

-2) 

0.805 

No 

2.80 

190 

0.366 

0.117 

6.5 

0.012 

51 

369 

90.0 

6.5(- 

-2) 

1.389 

Yes 

0.76 

10.4 

2.503 

0.146 

2.5 

0.008 

27 

302 

5.42 

1.1(- 

-1) 

0.664 

No 

4.17 

392 

0.307 

0.109 

5.0 

0.008 

57 

366 

80.8 

1.7(- 

-1) 

1.795 

Yes 

0.93 

7.48 

5.937 

0.128 

6.0 

0.008 

68 

374 

89.6 

l-2(- 

-1) 

1.197 

Yes 

1.40 

8.80 

5.799 

0.135 

5.0 

0.004 

81 

377 

84.4 

2.2(- 

-1) 

1.141 

Yes 

3.62 

10.8 

16.79 

0.138 


Computed with Reimers mass loss on the AGB. All other models have Vassiliadis & Wood (1993) mass loss. 
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Table 2: Surface abundance results for Zn through to Sr, taken at the tip of the AGB from 
the last computed model. 


Mass PMZ [Zn/Fe] [Ge/Fe] [Se/Fe] [Br/Fe] [Kr/Fe] 5 86 Kr/ 82 Kr a [Sr/Fe] 

A Fe b 

A = 0.02, Anders & Grevesse (1989) C, N and 0 



3.0 

0.0 

0.001 

0.004 

0.004 

< -0.001 

0.006 

-56.05 

0.007 

< 0.001 

3.0 

0.001 

0.071 

0.263 

0.346 

0.299 

0.477 

-693.3 

0.700 

-0.00152 

3.0 

0.002 

0.122 

0.397 

0.493 

0.433 

0.628 

-774.4 

0.853 

-0.0021 

4.0 

0.0 

0.004 

0.011 

0.009 

0.003 

0.010 

-25.82 

0.007 

-0.0016 

4.0 

1(— 4) 

0.077 

0.282 

0.365 

0.282 

0.519 

-40.31 

0.587 

-0.0026 

5.0 

0.0 

0.012 

0.031 

0.026 

0.011 

0.028 

-10.85 

0.011 

-0.0024 

5.0 

1(— 4) 

0.024 

0.096 

0.132 

0.083 

0.241 

307.7 

0.246 

-0.0024 

5.0 C 

0.0 

0.073 

0.188 

0.181 

0.107 

0.210 

91.56 

0.085 

-0.0090 

6.0 

0.0 

0.038 

0.106 

0.104 

0.051 

0.131 

92.09 

0.055 

-0.0038 

6.5 

0.0 

0.047 

0.133 

0.138 

0.064 

0.187 

186.3 

0.083 

-0.0037 




A = 

0.012, Asplund et al. (2005) C, N and O 




3.0 

0.002 

0.106 

0.434 

0.566 0.473 

0.645 

-747.2 

1.232 

-0.0028 

6.5 

0.0 

0.042 

0.113 

0.111 0.054 

0.136 

55.60 

0.053 

-0.0047 

6.5 

1(— 4) 

0.052 

0.164 

0.203 0.122 

0.397 

1171 

0.300 

-0.0048 

Z = 0.008 

2.5 

0.002 

0.148 

0.478 

0.600 0.547 

0.794 

-722.8 

1.31 

-0.0040 

5.0 

0.0 

0.163 

0.380 

0.401 0.302 

0.506 

355.3 

0.264 

-0.0138 

6.0 

0.0 

0.179 

0.432 

0.476 0.324 

0.657 

712.7 

0.388 

-0.0098 

Z = 0.004 

5.0 

0.0 

0.403 

0.813 

0.897 0.734 

1.176 

1178 

0.927 

-0.0216 


a 5( l X/' cei X) notation is defined according to \{ l X/' [ei X)/{ l X/ ret X)(T ) — 1] x 10 3 where ref A is a reference 
isotope. 

^Defined by log(lfi n ai/Yinitiaih where Y is the surface abundance in number fraction. 
c The model has Reimers mass loss on the AGB instead of Vassiliadis & Wood (1993). 
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Tabic 3: Abundance results for 0, Ne, P, S, Cl. and Ar. 


Mass 

Z e(0) a e(Ne) 

Ne/O b e(P) 

e(S) e(Cl) e(Ar) 

5 38 Ar/ 36 Ar c 

initial 

0.02 8.935 8.103 

0.181 5.591 

7.277 5.143 6.570 

0.000 


3.0 

no PMZ 

8.940 

8.279 

0.273 

5.595 

7.299 

5.192 

6.592 

7.735 

4.0 

no PMZ 

8.933 

8.204 

0.234 

5.613 

7.294 

5.190 

6.586 

14.55 

5.0 

no PMZ 

8.931 

8.191 

0.227 

5.617 

7.310 

5.208 

6.603 

17.76 

6.0 

no PMZ 

8.908 

8.168 

0.227 

5.696 

7.328 

5.237 

6.621 

26.29 

6.5 

no PMZ 

8.905 

8.166 

0.228 

5.635 

7.330 

5.212 

6.623 

25.18 


initial 

0.012 

8.659 

7.838 

0.189 

5.358 

7.139 

5.498 

6.178 

0.000 

3.0 

no PMZ 

8.656 

8.086 

0.336 

5.393 

7.165 

5.530 

6.205 

17.00 


0.002 

8.666 

8.666 

0.486 

5.522 

7.165 

5.534 

6.204 

41.34 

6.5 

no PMZ 

8.563 

7.868 

0.277 

5.410 

7.188 

5.548 

6.227 

61.39 

initial 

0.008 

8.514 

7.675 

0.181 

5.155 

6.849 

4.715 

6.142 

0.000 


2.5 

no PMZ 

8.521 

7.846 

0.431 

5.178 

6.870 

4.754 

6.162 

16.37 


0.002 

8.528 

7.985 

0.653 

5.251 

6.867 

4.762 

6.160 

32.81 

5.0 

no PMZ 

8.419 

7.817 

0.313 

5.454 

6.899 

4.886 

6.192 

92.67 

6.0 

no PMZ 

8.317 

7.781 

0.364 

5.471 

6.908 

4.861 

6.201 

63.44 


initial 

0.004 

8.206 

7.367 

0.181 

4.847 

6.542 

4.407 

5.835 

0.000 

5.0 

no PMZ 

8.052 

7.665 

0.512 

5.444 

6.609 

4.633 

5.900 

136.7 


a e(Y) = log 10(Y/ if) + 12, where abundances are measured by number, and H is the abundance of hydrogen. 
h The Ne/O ratio is calculated by mass, to allow comparison to the observations. 
c The 5 notation was defined previously in Table 2. 
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Bipolar PNe 


Bipolar PNe 



Fig. 1. — Histograms of Se and Kr abundances from SD08, separated into 0.1 dex bins, for 
non-bipolar, non-Type I PNe (top); Type 1 PNe (middle); and bipolar PNe (bottom). Data 
are shown only for objects in which Se and/or Kr emission was detected. Type I and (to a 
lesser extent) bipolar PNe exhibit smaller enrichments than PNe with less massive progenitor 
stars. 
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Fig. 2. — The Se and Kr abundances of PNe plotted against He/H (top panels), N/0 (middle 
panels), and C/0 (bottom panels). The triangles correspond to Type I PNe, the boxes are 
non- Type I bipolar PNe, and circles represent all other PNe from the sample of SD08. Se and 
Kr upper limits are depicted as open symbols. Representative error bars to the abundances 
are shown in the left hand panels. The observed correlations show that Type I PNe generally 
display smaller s-process and C enrichments than other PNe. While many bipolar PNe follow 
the same trends, some are significantly enriched in C, Se, and Kr, suggesting that they arise 
from less massive stars. 
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Fig. 3. — The temporal evolution of the (a) hydrogen (solid line) and helium-exhausted 
(dashed line) cores, (b) the He-shcll temperature and (c) the temperature at the base of the 
convective envelope for the 6.5 M 0 , Z = 0.012 model during the TP-AGB. 
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Fig. 4. — We show the [Se/Fe] ratio (top panel) and the [Kr/Fe] ratio (bottom panel), as a 
function of the thermal pulse number, from a selected number of stellar models. Each point 
represents the surface abundance during the interpulse period. 
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Proton number, Z 

Fig. 5. — The surface abundances of Zn, Ge, Se, Br, Kr, and Sr at the tip of the AGB phase 
for a selection of the stellar models. Abundances are shown as [X/Fe] ratios, and are plotted 
as a function of proton number, Z . The 3M 0 , Z = 0.02 model includes a PMZ of O.OO2M 0 
whereas all other models are without 13 C pockets. The gray-hatched box with a maximum 
at [X/Fe] = 0.3 represents the region of observed Type I PN Se and Kr abundances. Note 
that observational uncertainties allow for [Se, Kr/Fe] up to 0.6. Boxes around the elements 
Se and Kr indicate that for the present study we can only use these elements for comparison 
as large PN abundance data-sets of the other trans-iron elements are not yet available. 


